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Electron relaxation in quantum dots is studied theoretically in polar semiconductor materials, 
with an emphasis put on the phonon-bottleneck problem and the electron-LO-phonon coupling. The 
theory is based on multiphonon states of the electron-phonon system and the self-consistent Tamm- 
Dancoff approximation is used for the electronic self-energy. Electronic relaxation rate is shown 
numerically to be on the scale from hundreds fs to tens ps, for electron energy-level separations 
being in the broad range from about one LO-phonon energy to about three or four optical-phonon 
I energies. Despite of displaying some resonance features, the electronic relaxation rate does not 

. appear to be crucially dependent on the quantum dot size. 
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I. INTRODUCTION 



The semiconductor quantum dots (QD) have attracted attention due to the earlier expectations 1^,^ concerning 
■ their use in high efficiency semiconductor lasers. These expectations of the high efficiency were based mainly on 
the (5— function like electronic density of states of such zero-dimensional (OD) structures and on an expected narrow 
' gain region. In later experiments the luminescence of these OD structures was found to decrease with decreasing the 
^— j , lateral size of quantum dots |^-||] . Besides attributing this luminescence decrease to technological difficulties in the 
process of realization of the lateral confinement for electrons and holes, a different explanation was suggested 
ascribing the luminescence decrease to the so called 'phonon-bottleneck effect' (see ref. 0). The phonon-bottleneck 
^ hypothesis was supported by an earlier detailed analysis of the electron-phonon interaction in quantum dots based on 

d ' Born approximation to electron-phonon scattering ]q-pd[| . 

^ I In experiments, the relaxation of electronic energy was often measured on the scale of picoseconds fl^-^, being 
I ' therefore fast enough to express doubts about the existence of the bottleneck effect. A significant dependence of the 
^ , electronic relaxation rate on the quantum-dot size was not reported in papers |p^-^ and the relaxation efficiency 
g ' was reported to be rather independent of the relation between the electronic energy-level separations and the optical- 
^ ] phonon energy |l^ . Basing on the experimental observations, it has been pointed out that the relaxation process of 
electrons in quantum dots should be considered as a multiphonon process [p| p^|p^Jl^ -pl[| , and that the multiphonon 
transitions should be expected to be the main mechanism providing the electron-energy relaxation. 

The quantum dot system appears to be relatively simple because the electronic structure of the bound states, 
' unperturbed by the lattice motion, may consist of only several discrete energy levels. It is well known that the 
electron-phonon system of the quantum dot can be exactly diagonalized when the coupling of the electrons to the 
lattice vibrations is restricted to the transverse coupling (see below) terms only Another simple and exactly 

soluble model, with no dissipation of the electronic energy, is the system with two electronic energy levels coupled 
to a single mode of the lattice vibrations , with the Hamiltonian formally equivalent to Jaynes Cummings model 
A similar non-dissipating electron-phonon system was studied in a one-dimensional system |^5| . Although we 
are aware of several simplified systems in which the electronic subsystem does not relax the energy, it is not yet clear 
which of their properties can be generalized to a little more complicated realistic systems which we likely meet in the 
real quantum dots, namely in the systems, in which the carriers, confined in the dots, interact with a large number 
of the lattice vibrational modes, and in which the electron-phonon coupling is more general than that of the Jaynes 
Cummings model. Unless a further progress is made in the area of solving the electronic relaxation problem in the 
zero-dimensional structures by exact methods, the approximative theoretical approaches are in order. 

The role of the multiphonon states in the formulation of the electron transport problem in quantum dots was 
emphasized in papers |7[|2^-|30[| . The need for the self-consistent treatment of the effect of the " collision broadening" 
of the electronic energies in quantum dots was emphasized in references |2^j2^ . These two requirements were recently 
taken into account in the study of the electronic spectral density ||7|,^9|. Very sharp spectral density peaks were 
obtained. This sharpness was interpreted as an indicator of a very long electronic lifetime, without paying a sufficient 
attention to the lineshape of the spectral density features. In our previous paper pit] it was shown that the spectral 
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density peaks may diverge to infinity as the inverse of the square root of the energy variable. Although such maxima 
in the electronic spectral density function are integrable, the relation of the "width" of such spectral features to the 
electronic lifetime may not be as simple as it may be in the case of the Lorentzian peak shape. The relation of 
the electronic spectral densities to the rate of relaxation of the electronic energy, reported preliminarily in references 
is treated in the present paper. 

Besides the electron-LO-phonon coupling, other mechanisms were considered recently as possible candidates for the 
explanation of the experimental data on the fast electronic relaxation in quantum dots. So, the finite lifetime of the 
optical phonon was shown to provide an efficient mechanism of relaxing the severe restrictions imposed by the energy 
conservation in the Born approximation upon the electron-LO-phonon scattering [ p4|j3^ . Also, the ultrafast electron 
energy relaxation in quantum dots has been recently suggested to be explained by the interaction of the carriers with 
the defect states in the quantum dot, taking into account the lattice relaxation mechanism | p6[ |. 

It has been shown recently, that the electronic scattering in the low-dimensional structures should be treated 
with caution. Namely, the Born approximation to the electronic scattering, giving good results in the Monte Carlo 
semiclassical simulations of the electronic transport properties of the bulk semiconductor samples, appears sometimes 
to be rather insufficient in the quasi-two and quasi-one dimensional structures like quantum wells and quantum wires 
p?! , ^ . These observations provide additional arguments in favor of going beyond the Born approximation in the 
electron-phonon scattering in quantum dots. Understanding the failure of the Born approximation as an implication of 
the multiple reflections of the charge carrier from the boundaries of the low-dimensional structure, it can be expected 
that such dimensionality effects should be present in the zero-dimensional structures like quantum dots. In such a case, 
the multiple scattering of electrons should be taken into account and the multiphonon states of the electron-phonon 
system of the quantum dot should be of importance. 

The purpose of this paper is to develop further the theory of the electronic relaxation in quantum dots, basing 
on the multiple phonon scattering of the electrons, considering the interaction of the quantum dot electrons with 
the longitudinal optical phonons of the bulk matrix sample into which the quantum dot is built, and to present the 
numerical results of the electronic relaxation rate in dependence on the separation of the electronic energies and on 
the temperature of the lattice. 



II. THE MODEL AND HAMILTONIAN 



Generally, there may be several electrons and several holes in a quantum dot structure. Assuming the charge 
neutrality of the system one can speak then about single-exciton, bi-exciton , etc., states, which are strongly influenced 
by the carrier-carrier coupling, expected to be signiflcant in quantum dots |39||40(| . The inclusion of the carrier-carrier 
coupling may complicate the analysis of the carrier relaxation in quantum dots, although it may be rather significant 
from the point of view of a quantitative comparison with experimental data. It is the purpose of this paper to 
concentrate the attention on the electron-phonon mechanism of the carrier relaxation. Therefore, the simplest model 
system is chosen to be considered, namely that consisting of a single electron in the quantum dot, while the presence 
of holes in the valence-band states is completely ignored. The authors believe that the basic features of the fast 
relaxation mechanism of the excited carriers in quantum dots, based on the electron-LO-phonon Frohlich interaction, 
are contained already in this relatively simple single-electron model. The inclusion of the carrier-carrier correlation is 
therefore left to a future work. 

In order to simplify the numerical part of the work, the quantum well is assumed to be of cubic shape, with infinitly 
high electronic potential energy outside the well. The penetration of the electronic stationary wave functions into 
the potential barriers is therefore neglected. Conflning the electronic motion to the lowest lying electronic states is 
believed to minimize the effect of this infinitly-deep quantum-well approximation. The Schrodinger equation with the 
electronic effective mass inside the well gives a set of solutions. From these electronic states we consider only two 
states in the present model: the ground state ipo, with the energy Eq, which will be put equal zero here, and one 
excited state V'l belonging to the triply degenerate first excited energy Ei. This set of two nondegenerate electronic 
unperturbed eigenstates provides a minimum basis set for considering the electron-energy relaxation due to electronic 
transitions between two electronic states coupled via the electron-phonon interaction. The complete electron-phonon 
Hamiltonian H then reads ||3l[] 

H = Ho + Hi, (1) 

where 

1 

Ho^Yl Enc+Cn + J2 ^Lob+b^, (2) 

n— q 
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and Hi is Frohlich coupling 



m,n—Q q 

of the electron to the system of the dispersionless optical phonons of the whole bulk of the sample with volume V, 
inside of which the quantum dot is built. In the latter equation q =\ q\, Aq = (— ie/g)[£'i,o('«^^ — /^o ^)]^^^(2£ol^)~^^^, 
where Koc and kq are high-frequency and static dielectric constants, Eq is permittivity of free space, e > is electronic 
charge. The operator bq annihilates the LO phonon with wavevector q, while c„ annihilates electron in state ipn- The 
electron-phonon coupling depends on the quantum dot size via the form-factor 

$(n,m,q)= J d'rrn{r)e''"'^^{r), (4) 

where the integration spreads over the quantum dot volume. It has been demonstrated earlier that in the case of 
calculating integral quantities like the electronic relaxation rate, it is rather plausible to neglect the impact, which the 
interfaces in a quasi-two dimensional heterostructure of GaAs-Ala;Gai_2: As-type might have on the calculated results 
^2| . It is assumed that the quantum dot structure is based on the GaAs-A^Gai-^: As-type heterostructure and that 
the lattice dynamics of this structure can be approximated well enough by that of the bulk GaAs. In other cases, like 
in CdTe quantum dots dispersed in poly (vinyl butyral)(PVB]_H], the influence of the interface modes of the optical 
lattice vibrations may deserve a more detailed attention ||4^ , |45| . 

Let us remark, that the operator Hi of the electron-phonon coupling contains two kinds of terms. The terms with 
n = m express such an interaction of the electron in the state n, in which the electron emits or absorbs a phonon, but 
the state of the electron remains unchanged. These term are called transverse interaction terms. The longitudinal 
interaction terms, with n m, are those which lead to a change of the electronic state upon absorbing or emitting a 
pho non ]4^ . When the longitudinal terms are omitted, then the remaining Hamiltonian can be diagonalized exactly 
|p2| . Although the longitudinal terms may appear to be the only important ones for the electronic scattering process 
of the electronic transfer from the excited state to the lower energy state, with emitting or absorbing simultaneously 
the phonon, both the longitudinal and the transverse terms will be taken into account in this work. The interaction 
operator Hi considered in this work does not have the form of the Jaynes Cummings operator and the electronic 
system is coupled to a large number of the phonon modes, so that approximative theoretical methods will be used in 
the following analysis. 



III. RELAXATION RATE 



In photoluminescence experiments an electron may be excited by a light pulse from the valence-band states into the 
conduction band states and can be finally prepared in an electronic excited state in the quantum dot. In the case of 
a time-resolved detection of the luminescence, corresponding to the process of annihilation of the electron with a hole 
in the valence-band states, with the simultaneous emission of a light photon, the rise time of the luminescence can be 
determined experimentally [p^-p^ . The value of the rise time is limited from below by the time which is needed for 
the transport process in which the electron is transferred from the excited state to the electronic ground state in the 
dot, from which the electron recombines with a hole in the valence band states. In terms of the present simple model 
of the quantum dot one may say, that in this way the time-derivative dNi/dt of the population Ni =< c^Ci > of the 
electronic state V'l can be related to the experimental data. 

In the present work, the latter transfer mechanism of the electron, between the electronic excited and the ground 
state, is supposed to be the process of the scattering of the electron on the system of the optical phonons. In the 
model of the quantum dot, considered here, it is assumed that at time t = the system is prepared in the state 
with the single electron occupying the excited state only, while the phonon subsystem is found at equilibrium with 
temperature T^. 

The process of the electronic relaxation can be theoretically formulated within the theory of the nonequilibrium 
statistical operator |^^. The rapidity of the relaxation of the electronic energy can be expressed with the help of 
the time derivative of the population, d < ci > /dt, of the electronic excited state {n ^ 1). In the nonequilibrium 
statistical operator theory, assuming not too short time scale, the state of the system is assumed to be described by a 
set of quantities < Pk >, with the corresponding set {Pi, P2, ■ ■ ■ , Pk, ■ ■ ■} of operators, suitable for the description of 
the system. The averaging in < Pfc > is performed with the nonequilibrium statistical operator. In the case presently 
considered the electronic subsystem will be characterized by the mean values Ni =< c^ci > and Nq =< CqCq > 
giving the average population of the electronic excited state and the ground state, respectively. The state of the 
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phonon system, assumed to be permanently at the thermodynamic equilibrium, will be characterized by the mean 
value < 6q6q > given by the Bose-Einstein distribution function at the lattice temperature T^, independent of time. 
Confining our choice of {Pm} to < c^ci >, < CqCq > and < foqfoq >, the kinetic period of the evolution of the 
nonequilibrium system under study is presumed. 

The "well-known expansion of the nonequilibrium statistical operator in powers of the interaction leads to the well- 
known expansion of the " collision integral" in the generalized kinetic equation giving the time evolution of the mean 
value < Pfe > . The lowest-order (in Hi ) terms of the collision term of this equation then read : 

where 

sj^^ = ^<[P,,Ho]>l (6) 

and 

sl^^^^<[P,,Hi]>l. (7) 

f(2) 

The term S'j, , which is at least of the second order in Hi , contains generally a factor expressing the memory properties 
of the system | |4^ ] . In the markovian approximation p7[ | , considered here, it simplifies to the following expression: 



= f dtie^'' < [Hi{ti), [Hi.Pk]] >l . (8) 

1^ J — oo 



The averaging < ... >q= Tr{. . . pq{t)) is performed with the so-called quasi-equilibrium statistical operator pq 
defined as 



Pi 



it)^Q-'exp{-J2Fk{t)Pk}. (9) 



Qq^Trexp{-Y,Fk{t)Pk}, (10) 

k 

in which the quantities Fk(t), playing the role of the intensive quantites describing the system at equilibrium (e. g. 
temperature), are determined from the conditions 

< Pk >l=< Pk > ■ (11) 

Consistently with our choice of the operator set {Pi, P2, Pk, ■■■}, the term s\^^ is zero. The coupling of the 
electronic subsystem to the laser light, which is not considered here, may lead to the appearence of the interband 
electronic polarization, which would be then a manifestation of the coherence between the electronic and light systems 
p8[ . Similarly, the field of the LO phonons can be coherently coupled to the system of excitations of the electronic 
subsystem. These effects could be obtained already from the term Sj^K The right-hand side of eq. (7) gives then the 
average value of the product of two electronic and one phonon operators. Consistently with our choice of the operator 
set {Pk} the term S*^^-* will be ignorred. Extending however the operator set in a suitable way, the Rabi oscillations, 
studied recently by Inoshita et al. |^^, could be considered simultaneously with the term S^^K The coherent coupling 
between the electrons and the LO phonons may play a role especially in the range of short time after the light pulse. 
This effect will not be considered in this work. The lowest-order part of the time derivative of < Pk > is then given 
by the formula (||). 

In the formula (|^), giving the lowest-order markovian contribution to the relaxation rate d < Pk > /dt, the 
operator Hi(ti) is the interaction operator Hi expressed in the Heisenberg representation. When approximating this 
Heisenberg representation by the interaction representation, the formula (^) leads directly to the well-known Born 
approximation to the collision integral. On the other hand, the Heisenberg representation of the interaction operator 
Hi in (||) makes it possible to sum partially the terms contributing to S'^"^^ up to infinite orders in Hi | |47[ |. 

The emphasis will be put now on the obtaining the kinetic equation for the time evolution of the population 
< c^ci > of the upper electronic state. At first, the commutator [i?i,c^ci] is calculated. As a result, several terms 
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are obtained, each being a product of two electronic and one phonon operators. The average in the formula then 
reads: 

< [i7i(ti),[i?i,c+ci]] A,Ap^r,s,p) (12) 

r,s,n,q,p 

X { $(n, 1, q) < [(bpih) - btpih))c+{h)cs{h), (&q - 6tq)c+ci] >l 
-$(l,n,q)< [(6p(ii)-fetp(<i))c+(ti)c,(ti),(feq-6tq)c+c„] >* }. 

Writing down the commutators in the latter equation explicitly according to the definition, the average quantity 
A =< [Hi{ti), [Hi,c^ci]] >* appears to consist of a number of terms. Each of these terms is an average of the 
product of six particle operators, three of them being in the Heisenberg representation. An example of these terms 
contributing to A is: 

J2 ApA^ir, s, p)$(l, n, q) < 6tp(ti)c+(ti)c,(ii)feqC+c„ >^ . (13) 

r,s,n,p,q 

The average value in the latter formula, and in the other terms contributing to A, will be decoupled according to the 
following scheme: 

< btp{ti)bqc+{h)cs{ti)ctcm >* ~< foip(ii)&q >*< c+(ti)c„ >*< Cs{ti)c+ >\, (14) 

in which only the particle operators belonging to different interaction operators are paired. All the anomalous averages 
of the type hke < 6p (ti)6q > and < c+(ti)c+ > are considered to be zero. The latter kind of the decoupling of the 
averages of the operator products in the relaxation rate formula was studied earlier in connection with the electron 
energy collision broadening and the electronic relaxation [ ^9|j5C| ] in bulk samples. It was shown, that the above type 
of the decoupling leads to the relaxation rate formula with plausible thermodynamical properties. Although this 
decoupling scheme appears to work well in the bulk structures, it is not obvious whether it is suitable enough in 
the presently considered two-level zero-dimensional system. This decoupling is only an approximation to the full 
correlation function (13). The question of the validity of this decoupling will not be studied in this work. 
In the above obtained single-particle correlation functions, like < Cs{ti)c^ >*, the diagonal approximation 

< Cs{tl)ct >\ ~ < Cs{tl)4 >\ Ss,r (15) 

is assumed to hold. This assumption means that the quantum number n, indexing the unperturbed electronic states 
il^m remains a good quantum number and that the effects of the quantum coherence, included eventually in terms like 
< ci{ti)cQ >*, do not change considerably, in materials with a rather weak electron-phonon coupling like GaAs, the 
resulting picture provided by the diagonal approximation. The way of decoupling the average A is consistent with 
the above-made assumption about the choice of the operators Pk describing the system under study. Nevertheless, 
the possible influence of the nondiagonal terms like < ci{ti)cQ >^ remains to be verified in a future work. 

The phonon correlation functions will be taken into consideration in the diagonal approximation too. This approx- 
imation would not be appropriate when considering the hot phonon effect with the LO phonons produced in the area 
of the quantum dot. The hot-phonon effect will be neglected here. Assuming that the phonon one-particle correlation 
functions like < 6q6q(ti) > are invariant under the operation of q — s- — q, we get: 

<[H,{h),[H^,c+c{\]>\ (16) 

^2ReY,\ A, n $(1,0, q) p { (< 6q(ii)&+ >* + < 6+(ti)6q >*) 
q 

X(< C+{ti)ci >*< Co(tl)4 >g - < Ci{ti)4 >\< 4{tl)cQ >*) }, 

where Re denotes the real part of a complex number. The knowledge of the time ii-dependence of the single-particle 
correlation function in the latter equation would make it possible to perform the integration over ti in (||) and to get 
in this way the rate of change of the population < c^ci >* as a function of the time t. 

In calculating the single-particle correlation function we shall proceed in the following way: The single-particle 
correlation function will be determined formally for the system with the Hamiltonian (1), in the case of the ther- 
modynamic equilibrium at a temperature T and the obtained functional dependence of the electronic and phonon 
quantities on the temperature will be formally transferred to the case of electrons and phonons having different 
respective temperatures and T^. 
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In the fo llowin g the electronic correlation functions will be expressed in terms of spectral densities (see eqs. ( A12| ), 



( |A13| ) and ( A14 ) in Appendix A). The electrons and phonons will be assumed to have different temperatures. Namely, 



at time t — Q with A^i = 1, the phonon distribution function ly^o —< &q^q > will be given by the Bose-Einstein 
distribution function taken at a temperature T^. In the two- level electronic system, unperturbed by the electron- 
phonon interaction, with the energy levels Eq and Ei, the chemical potential is fj, = {Eq + Ei)/2. The population A^i 
of the excited electronic state at a temperature Tg is 

^1 = 7p-^' (17) 

(ks is Boltzmann constant). For example, the state of the electronic subsystem with one electron in the excited 
state and with the empty ground state then corresponds to the limit of Tg ^ 0_ . In this way the population of this 
two- level system can be formally expressed with the help of the electronic temperature Te. For the purpose of the 
present calculation it is suitable to describe the state of the two-level electronic system in terms of the population of 
the electronic state rather than in terms of an electronic temperature. 

Considering the process of the electron-energy relaxation as a markovian process and taking into account the 
decoupling (16), the theory of the nonequilibrium statistical operator leads finally to the following formula for the 
rate of change dNi/dt of the population of the electronic state with n = 1: 



dN lir ' { 

^ =-— am Ni{l- No) { {I + VLo) dE ai{E)cjo{E - Elo) (18) 



dt h 



/oc 
dEai{E)(Jo{E + Elo) 
-OO 

-No{l - Ni) ( (1 + lyLo) I dEao{E)ai{E - Elo) 

\ J — OO 

/OO 
dEao{E)a^{E + ELo) 
-oo 



where ctq and cti are electronic spectral densities (see eq. ( A14)), v^o is Bose-Einstcin distribution of LO-phonons at 



temperature Tl of the lattice. The constant amn is defined as 

ttmn = ^ Mg Pi $(?l,™,q) P, a-mn = anm- (19) 

q 

The first two terms in the square brackets in eq. ( p^ ) are the only terms contributing at the the initial moment 
of time {t = 0), at which Nq = 0. The first term provides the production of phonons, while the second term has 
the meaning of decreasing the electronic population Ni while absorbing an LO phonon. This second term, with the 
absorption of a phonon, appears in the formula from two reasons: First, the relaxation rate formula (18) does not 
contain the energy conservation (5-function as it would appear in the first order of the time dependent perturbation 
calculation formula. Second, the term with the absorption of a phonon depends on the overlap of the two spectral 
densities, cti and do- This overlap of ui and do is determined by the properties of the interacting electron-phonon 
system, namely, by the multiplc-phonon nature of the eigenstates of this system, which may allow for the non-zero 
contribution of this phonon- absorption term to the relaxation rate at t = 0. In fact, the presently used approximation 
to the electronic self-energy (see below) assumes the multiple-phonon property of the eigenstates of the system. The 
numerical values of the contribution of the phonon-absorption term comes out nevertheless rather small in GaAs (see 
below). This is in agreement with the fact, that the electron-LO-phonon coupling in the bulk GaAs is not strong. 
The last two terms in the square brackets have the meaning of a transfer of the electron from the state n = to the 
excited state n = 1. These terms would play a role at the later stages of the relaxation, namely at i > 0, which is 
however not treated here. 

The electronic self-energy is determined in the self-consistent Tamm-Dancoff (or the self-consistent Born) approx- 
imation [Q. The use of the self-consistent Tamm-Dancoff approximation makes it possible to include the states in 
which the electron is coupled various numbers of phonons and to go beyond the Born approximation when determin- 
ing the self-energy of this highly singular system with (5— function type unperturbed electronic density of states. It is 
shown in the Appendix B that after a simplification, concerning the electronic distribution function Nm, the Dyson 
equation of the (retarded) electronic self-energy for the electron in the state n can be brought to the following form 
ill: 
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1 

(20) 

m— 

^ r 1 - jV^ + VLO 

E-E^ + Elo- M„,{E + Elo) + *0+ ^' 

where iV^ is the electronic population of the m-th state. In the numerical evaluations below it is assumed that the 
system is at the initial state with TVi = 1 and Nq = 0. The above Dyson equation for M„(£') assumes that the 
electronic Green's function, being generally a matrix in the index n in the representation of the states ipm can be 
approximated upon taking the nondiagonal terms as zero. This assumption means that we regard the index n to be a 
"good quantum number". Also this means, that we assume that the mixing of the electronic states V'l and V'o via the 
operator Hi does not contribute seriously to the leading terms of the resulting relaxation rate. This approximation to 
the electronic Green's functions is in agreement with what was assumed above concerning the correlation functions. 
The extent of the validity of such approximations in a quantum dot structure can be decided by a more detailed study. 



The spectral densities (T„(-E) can be obtained (see eq. (A14)) from the real and imaginary part of the retarded 
self-energy M„(£'), M„(£') = Rn{E) — iYn{E), Yn > 0, which is determined by the equation eq. (|2^). The derivation 
of this equation, together with the approximations applied, is presented in Appendix B. The solution of equation ( [20| ) 
can be performed in two ways: (i) after some simplification this equation allows for an analytical solution, (ii) the full 



numerical solution of eq. (20) can be performed in the present model. 



IV. AN ANALYTICAL SOLUTION 

An approximation can be made, under which the Dyson equation (|2^) can be solved analytically. This approximation 
consists of several simplifications, as follows: 

The first approximation step is the neglection of the real part of the self energy. The relaxation rate ( p^ is 
proportional to the product of the spectral densities, which depend significantly on the imaginary part of the self- 
energy. It will be seen later, that the direct comparison with the result of the numerical solution of the full equation 



( pop shows, that neglecting the real part of self-energy does conserve important properties of (20). 

On the basis of a simple Born approximation approach, it may be expected that for the relaxation process it holds 
approximately that an electron with energy Ei will make a transition to the state with the energy Eq and emits an 
optical phonon with the energy E^o and that the transverse part of the electron-LO-phonon coupling can therefore be 
neglected. In fact, the constants aoo and an characterizing the transverse coupling are larger than aoi characterizing 
the longitudinal interaction. The dependencies of amn on the detuning between the electronic excitation energy Ei 
and the optical phonon energy E^o are shown in Fig. |l] for GaAs. It may also be expected that the imaginary part 
of the self-energy, Yn{E), is significantly nonzero near the energy oi E = En- Therefore, one of the terms on the rhs 
of equation (|2^) gives the leading contribution, while the other term can be neglected. Under these approximations 
the Dyson equation ( po|) can be written in the form of the following set of equations for Yn{E) (reminding that 
MniE) ^ Rn{E) ~ lYniE)): 

^ ' \E-Eo-ELo+iYoiE-ELo) + 0+}' ^ ' 

^ ' \E-Ei+ELo+iYiiE + ELo)+0+} ^ ' 

The real part of the self-energy was neglected in this equation. The non-negative solutions of this equation set are: 



^0 ^ /I ^ ' ' ^ ^^-\E\\E-Ei+Elo I, (23) 
\/ \ E - Ei+ Elo I 



Yi{E)^^I^^^^J^-\E-Ei\\E-Elo\, (24) 
\/ \ E - Elo I 
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where ^ = aoi(l + i^lo) and we have put Eq equal to zero. It is seen, that the imaginary part of the self-energy Yn 
goes to zero as y/E at the energy E — En- Also, Yi diver ges at = Elq and Yq diverges at E = Ei — Elq, both 
like I/^/E. The corresponding spectral densities (see ( A14 )) then diverge at these points on the energy axis like \/E. 

The plot of Yq{E) and Yi{E) was presented in the Rapid Communication ||3l[] and will not be repeated here. It is 
demonstrated there, that some peaks in (T„(£') of a state with a given n can be understood as being phonon satellites 
of the electronic state with m ^ n. This property of the self-energy and of the spectral density guarranties, that the 
two spectral densities, as they occur in the formula ( ]l8| ) for the relaxation rate, have always a nonzero overlap, no 
matter how large the quantum dot is, or, in other words, how large the detuning between Ei — Eq and E^o is- This 
property of the self-energy can be seen as a formal expression of the idea of the absence of the phonon bottleneck, 
even in the present case of electrons interacting with dispersionless phonons. 

The relaxation rate dNi/dt (at t > Q) computed in the approximation of the equations (22) and (23) to the electronic 
self-energy, is shown in Fig. |[ In this Figure the relaxation rate reaches the order of 1/ps in a range of the dot size, 
in which Ej^o — Ei varies in the range of about 30 meV. The relaxation rate in Fig. ^ demonstrates, that the present 
mechanism of the electronic energy relaxation, with the transverse electron-LO-phonon interaction neglected, can give 
the relaxation times in the range of picoseconds or hundreds of femtoseconds in a rather broad range of the detuning. 
The double-maximum shape of the relaxation rate in Fig. ^ is probably due to the approximations introduced in the 
analytical solution and due to the satellite structure of the electronic spectral densities obtained in this approximation 
to the self-energy Isij . 



V. NUMERICAL SOLUTION 



In this section the numerical solution of the full equation (20) is presented. That part of the electron- phonon 
coupling, which is proportional to the constant aoi, was shown in the previous section to give the effect of the absence 
of the phonon bottleneck and the relaxation times in the range of picoseconds. In the present section the coupling 
constants ago and an, characterizing the transverse coupling of electrons to the phonons, are taken into account and 
both the real and imaginary parts of the self-energy are included into the computation. 

The numerical calculations are performed for GaAs, InAs and InP |5^,Q. The Dyson equation ( pO| ) for the 
retarded self-energy is solved by iterations. In contrast to the previous section, the Eqs. ( pO| ) are solved numerically, 
substituting a positive finite number A instead of 0+ in (pO|). We used A = 10"'^ meV. The real and imaginary parts 
of the self-energy are presented in Fig. || for GaAs for the detuning chosen to be E^o ^ E^o = ^8 meV at Tl = 77 K. 
This Figure shows the rather complicated satellite structure of the electron self-energy obtained in the Tamm-Dancoff 
approximation. In the curves, displayed in the Figure, the phonon satellites of both the electronic ground and excited 
states are observed. 

The spectral density of the electron in the excited state would he 5{E — Ei) in the noninteracting case (-ffi = 0). 
In the interacting electron-phonon system the spectral densities <Jo{E) and ai{E) have a rather complex structure 
displayed in Fig. ^. In this Figure the numerical data are computed for the parameters of GaAs at the detuning 
Elo — El = —8 meV at Tl = 77 K. Namely, the optical-phonon energy Elo = 36.2 meV and the value of the energy 
of the excited unperturbed electronic state is chosen to be 44.2 meV. Both Figs. ^ and ^ show multiple peaks. It is 
seen that one of the characteristic energy separations among the individual graphs is 36.2 meV. This can be interpreted 
as follows: an electronic state in the quantum dot with the electron-phonon interaction appears to consist of such 
components, which correspond to the two unperturbed electronic states coupled to various numbers of LO-phonons. 
For example, the peak near E = in Fig. ^a corresponds to the phonon-less n = line, while the peak near 36meV 
corresponds to a one-phonon satellite of the electronic ground state feature (n = 0). In the same Fig. ^a the feature 
near 8meV can be seen as a one-phonon satellite of the excited state (n = 1), having the energy decreased from Ei 
by Elo- Similar energy relations are observed in the Fig. The main peaks in Fig. ^ appear to be spread over 
an interval of several meV. The shape of the peaks was characterized in the previous paper ||31| , where it was shown, 
that the shape of the individual maxima in the pattern of the spectral density on the quantum dot described by the 
present model, do not have the form of a Lorentzian. Although the present results do not allow to compare with the 
known experimental data on luminescence, it may be said that the occurence of the phonon satellites in the spectral 
densities is in a qualitative agreement with the experimental data on luminescence, in which luminescence maxima 
separated by one optical-phonon energy use to be observed (see e. g. p^j20| , ^ ). 

A considerable attention has been recently paid in experiment to the width of the luminescence lines p^-|56t in 
connection with the coherent phonons in the nanocrystals and semiconductor quantum dots. The issue of the variation 
of the "width" of a spectral line, whether it decreases or increases, with changing such parameters like the quantum 
dot size or the detuning, might help to judge the relevance of a particular theoretical model of the quantum dot or 
a nanocrystal. As it is seen in Fig. ^ and as it was discussed using analytical arguments in the reference the 
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shape of a feature in the spectral density differs from a Lorentzian, so that it not straightforward to speak about 
the hnewidth, at least in the case of the present model based on the LO phonons. An attention should be given 
to generalizing the present model in order to include the acoustic phonons and the electron-hole interaction. This 
question will not be further analyzed in this work. 

A significant property of the spectral densities <Jo{E) and cri{E) is their mutual overlap. Namely, from Fig. ^ it 
is observed that despite of the 'nearly discrete' structure of these spectral densities, the function cto(-E ± E^o) has a 
nonzero overlap with ai{E), in a rather broad and continuous range of the detuning Elo — Ei. This property of the 
spectral densities guarranties a nonzero electronic relaxation rate in the broad range of the detuning E^o ~ Ei. 

The numerical results giving the rate of the excited state depopulation, dNi/dt, as a function of the detuning 
Elo — El, are displayed in Fig. ^ for GaAs. The relaxation rate is computed for the initial state with TVi = 1 and 
A^o — and with two chosen temperatures of the lattice. The most intensive relaxation occurs, when the optical- 
phonon energy is at resonance with the excited-state energy Ei. In the region of small dots {Elo — Ei < 0) the rate 
displays an overall decrease with decreasing the dot size, at both temperatures. At such a detuning, which corresponds 
to El equal to 2Elo, ^Elo and AElo [Elo — Ei equal to about — 36meV, — 72meV and — lOSmeV, respectively) 
further resonance maxima appear. Although the sharp peaks of the spectral densities, as they are displayed in Fig. U, 
broaden when the lattice temperature is increased, the overall values of the relaxation rate decrease with increasing 
temperature. This decrease of the electronic relaxation rate with increasing the temperature of the lattice can be 
explained taking into account the nature of the electronic states determined in the self-consistent Tamm-Dancoff 
approximation to the electronic self-energy. The electron is coherently coupled to a number of the LO-phonons. This 
property of this electron-phonon system is responsible for the effect of a non-zero relaxation even at the condition of 
nonzero detuning Elo — Ei. As the temperature of the lattice increases, this coherence weakens, which leads to a 
decrease of the relaxation rate. 

It may be interesting to see, how much the process accompanied by the absorption of the LO-phonon (the term 
proportional to iVi(l — Nq)vlo in eq. (18)) contributes to the total relaxation rate at < = 0. The contribution of 
the latter term is given by the dashed line in Fig. ||, while the contribution of the main term, the one proportional 
to A^i(l — A'o)(l -t- VLo), is given by the full line. We see, that in accord with the fact that the electron-LO-phonon 
coupling in GaAs is rather weak, the phonon-absorption term gives only a rather small contribution to the total rate 
dNi/dt. 

The relaxation rate reaches the maximum of the order of about 1 electron per 100 femtoseconds. This value has the 
order of magnitude of the experimental observations in some studied samples |18|, in which the electronic relaxation 
rate was characterized by the relaxation time of several hundreds of fs. The relaxation rate computed at the room 
temperature of the lattice, at which the relaxation rate is decreased with respect to the low temperature case, is still 
at the scale of tens of picoseconds. These theoretical results, obtained within the present simple model of the quantum 
dot, therefore contradict the assumption of an existence of the phonon bottleneck in the electron-energy relaxation 
in quantum dots in polar semiconductors, at least in the case of small dots, such, in which the detuning is near to 
zero or negative. Let us remind, that the present theory is formulated for small dots (negative detuning). For large 
dots {Elo — i?i > 0) the present model has to be extended to include the effect of more unperturbed electron energy 
levels. This generalization of the model of the quantum dot is not done in this work. 

The overall efficiency of the presently considered mechanism of the electron energy relaxation depends on the 
strength of the Frohlich coupling. This dependence is presented in Fig. 0, in which the relaxation rate is displayed 
for InAs, GaAs and InP at 77 K of the lattice temperature. 

Summing up, using a single-electron model of quantum dot with two electron energy levels and with the electron 
coupled to dispersionless bulk optical phonons, this work presents the calculation of the electron-energy relaxation 
in quantum dots in polar semiconductors. The relaxation time is found to be in the interval from hundreds of 
femtoseconds to tens of picoseconds in such quantum dots, in a broad range of the detuning, in which the energy of 
the electronic excited state is equals from about one LO-phonon energy to about three or four times the energy of the 
LO-phonon. This overall agreement of the present results with the theoretical data suggests, that the electron-LO- 
phonon interaction provides an effective multiphonon mechanism giving the femtosecond or picosecond electron-energy 
relaxation time in a rather broad range of the quantum dot size of small quantum dots. Such conclusions apply to 
the temperature range from low temperatures of the lattice to the room temperature. 



APPENDIX A: EQUILIBRIUM CORRELATION FUNCTIONS 



In this Appendix several formulas are reminded, connecting the correlation functions with the spectral density and 
the self-energy in an equilibrium system at nonzero temperature. In particular, an approximate formula expressing 
the time dependent correlation functions in terms of the self-energy, suitable for a generalization to the case of a 
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nonequilibrium system, is presented. 

In the thermodynamic equihbrium at temperature T the time-dependent single-particle correlation functions can 
be spectrally decomposed in the well-known way |^ according to the formulas: 

< B{t')A{t) >= — / JBA{io)e'^'^' -'^dtu (Al) 

and 

< A{t)B{t') >= — / JBA(w)e*^e'"(*'-*)d^, (A2) 

where JbaI'^) is spectral intensity of the correlation function < B{t')A{t) >. The spectral intensity can be obtained 
from the knowledge of the retarded thermodynamic Green's function 



« A{ty,Bit') ^9it - t') < [A{t), B{t% >, (A3) 
in 



where the symbol < ... > denotes the averaging performed over the grand-canonical ensemble, 
A{t)=exp{iTi,t/?i)Aexp{—i'Ht/1l), Ti. = H — fiN, ji is the chemical potential and N is operator of number of par- 
ticles. Also, [A,B].^ = AB — -qBA, in which rj = ±1 for Bose and Fermi particles, respectively. Defining the Fourier 
picture <C A; i? | w as 

< A{t)- Bit') — / e-*'^^*-* ^dio « A; B I w (A4) 

and when the spectral intensity is real, then there is the relation between the Green's function and the spectral 
intensity 



Im <t: A:B 



\ uj :>r= -j^ [e'^'^B^ ~ T]^ Jba{uj)- ( A5) 



one obtains 



with 



Assuming in the present work that the influence of the elect ron-phonon interaction on the characteristics of the 
phonons is only minor, the phonon correlation functions will be determined in the zero-order approximation. Realizing 
that 

« 6q; b+ I uj J———, (A6) 

^ hio - Elo + «0+ 

< b^{h)b+ > + < fe+(ti)6q >= {j^LO + l)e-^'^-°'i + i^Loe'"^°*\ (A7) 

J^LO = . (A8) 

The electronic time-dependent correlation functions will be calculated with nonzero electron-phonon coupling Hi. 
In this case the spectral intensity will be expressed in terms of the electronic self-energy. In the diagonal approximation 
neglecting the inter-level correlations due to the electron-phonon coupling, the retarded Green's function of the n-th 
electronic state can be written as 

« c„; c+ I LJ »,= — ^ , (A9) 

where Mn \tiuj) is the retarded self-energy of the n-th electronic state. Working with the retarded self-energy only 

we shall drop the index (r) in the symbol Mn\hLu). Defining Rn{huj) = ReMn{huj) we shall write Mn{huj) — 
Rn(hLu) — iYn(huj). Such a decomposition of the self-energy into the real and imaginary parts allows one to work with 
the real and non-negative quantity Yn(huj) in what follows. 
With (|A2D the correlation function < c„(ti)c+ > reads: 
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< c„(ii)c+ >= - /" ) 

7^ J-oo \ e^B^ + 1/ 



— iujt\ 



duj. 



(AlO) 



Let us notice, that in the zero-order approximation, in which i?„ and Yn go t o zero, the expression expiUu) /ksT) in 
the latter formula can be substituted by exp((£'„ + n)/kBT). In formula (AlC) we shall approximately substitute the 
expression 1/ {exp(hu} / ksT) + 1) by the value of the Fermi-Dirac distribution function, giving the average number of 
electrons in the n-th electronic state. In this work therefore the electronic population will be approximated by 



(All) 



Two correlation functions then are: 



(A12) 



< c+{h)cn >= TiNn an{Lo)e"^''dw, (A13) 

J — OO 

where the spectral density (7n{E) is 

, ^ 1 K(M + o+ , . 

'^"^'^^ 7r(?iw-£;„ + Ai-i?„(?iL^))2 + (F„(M + 0+)2' ^ ' 

going to the (5-function in the hmit of zero self-energy and fulfiling the rule 

/OO 
CT„(tj)dw, (A15) 
-C30 

for n = 0, 1. 



The formulas (A7), (A12) and (A13) provide the dependence of the spectral densities on the temperatures, or on 
the populations, of the electrons and phonons, within the theory of the real time thermodynamic Green's functions. 
These functional dependences will be assumed to be valid for the correlation functions, which we meet in the formula 
( |r^ ) in the general case of nonequilibrium state of the system. In this way, the electronic relaxation rate (|l^) can be 
expressed in terms of the electronic spectral density and in terms of the electronic self-energy. 

The above given method of obtaining the single-particle correlation functions in the nonequilibrium system under 
study is accompanied by an inaccuracy, which may depend on the strength of the electron-phonon coupling. The 
reason for this is obvious upon comparing the statistical operators, with the help of which the averaging is performed. 
In the case of the therniodynamical equilibrium the averaging is performed in the grandcannonical ensemble of the 
system with the electron-phonon interaction included, while in the averages like < c^{t)cj >q, appearing in the 
formula (^^, the quasi-equilibrium statistical operator does not depend on the electron-phonon interaction, at least 
in the approximation we assumed in this work. This question can be clarrified with the use of a more systematic 
approach to this nonequilibrium system. 



APPENDIX B: ELECTRONIC SELF-ENERGY 

Using the language of the Feynman diagrams for the Matsubara Green's functions |5^j5^, the self-consistent Tamm- 
Dancoff approximation (or the self-consistent Born approximation) to the electronic self-energy can be expressed 
as a diagram containing two bare interaction vertexes connected with two lines, one of them corresponds to the 
full electronic Green's function G, while the other corresponds to the bare phonon Green's function D^^\ The 
equation for the self-energy of the electron in the n-th state, Mn{ihujp) , defined on the set of imaginary frequencies 
Up = {2n + 1)^!^^, n is integer, is 

MniiTiLjp) = kBTj2Yl I ^9 I'l *(™>".q) I' ^G™(«wp-iw,.)^^°n*^r). (Bl) 

m q r 
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Using the Lehmann's representation of the Green's functions |^ , the summation over the discrete imaginary frequen- 
cies LVr = 2r(jTkBT/h), r being an integer, can be performed. In this way the equation for the retarded electronic 
self-energy in the state m of is obtained: 

MniE) = J2 anrn[G:^'iE - ELo)il - npiE - Elo) + yw) (B2) 

+GZ\E + Ei,o){nF{E + Ei^o) + ^^lo)], 

where npiE) = 1 / {eyi\){E / K sTe) + 1) and vlq = I / {exp{ELo / ksTL) — 1) Here C*^* denotes the retarded electronic 
Green's function. Note, that in the latter equation the electronic Fermi-Dirac distribution function np{E) appears as a 
function of the energy variable E. Assuming the knowledge of the electronic temperature, this function would present 
no difhculty. However, in this work we approximate these electronic distributions by the Fermi-Dirac distribution 
function value Nm, which is the value of the Fermi-Dirac distribution function in the state with the unperturbed 



electronic energy Em, at which the corresponding retarded Green's function in the latter equation (B2) would have 
the pole in the case of no electron-phonon coupling. In this way the treatment of the electronic temperature can be 
avoided. 

The equation for the retarded self-energy is {Elq = TT'^lo)'- 



E-E,n~ Elo ~ M„,{E ~ Elo) + 

Nyn + ^LO 1 



E-E.^ + Elo - M^iE + Elo) + ^0^ 



(B3) 



The equation (|B3|) closely resembles the similar equation for the self-energy of Frenkel excitons studied earlier [ p9[ . 
In the reference |59(| this equation is derived with help of the real-time thermodynamic Green's functions. 

In (B3) the chemical potential is not written. This form of the equation corresponds to the choice of the Hamiltonian 
H in the equation of motion for the Green's function, instead of the Hamiltonian Ti. = H — fiN^ although the statistical 
averaging is still performed with the g rand -canonical ensemble. Consistently with this choice of Hamiltonian, the 



chemical potential fi is dropped in eq. (A_14) 
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FIG. 1. The dependence of the coupUng constants Uij in GaAs on the detuning between the optical phonon energy Elo and 
the excited state energy Ei . 
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FIG. 2. The dependence of the electronic relaxation rate {—dNi/dt) on the detuning Elo—Ei as computed for the analytical 
solution of the simplified Dyson equation (with the transverse coupling terms neglected) at the temperatures of the lattice 
Tl = 77 K. 
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FIG. 3. Real and imaginary parts of electronic self-energy in the states with n = 0, 1 for the detuning Elo —Ei = —8 meV in 
GaAs at 77 K. The full and dashed lines denote, respectively, the imaginary part of the self-energy (—ImMo) and (— /mMi), 
while the dashed (dotted) line denotes ReMo (ReMi). 
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FIG. 4. Spectral densities ao{E) (a) and cn{E) (b) in GaAs computed for the detuning Elo — Ei = — 8meV, at the lattice 
temperature of 77 K (full line) and at the room temperature (dashed line). Note that he E—axis of graph (b) is shifted by the 
energy of optical phonon with respect to graph (a). 
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FIG. 5. Relaxation rate {—dNi/dt) in GaAs as a function of the detuning Elo — Ei at two temperatures of the lattice. 
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FIG. 6. Two contributions to the total relaxation rate {—dNi/dt) in GaAs at Tl = 300 K. The full line gives the contribution 
of the process accompanied by the emission of LO phonon, while the dashed line gives the one of the process with LO phonon 
absorption. 
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FIG. 7. A comparison of the relaxation rate {—dNi/dt) at Tl = 77 K in three semiconducting materials. 
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